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Abstract 

We consider chiral perturbation theory with a nonzero 9 term. Because of the CP violating term, 
the vacuum of chiral fields is shifted to a nontrivial element on the SU{Nf) group manifold. The 
CP violation also provides mixing of different CP eigenstates, between scalar and pseudoscalar, 
or vector and axialvector operators. We investigate upto 0{6^) effects on the mesonic two-point 
correlators of chiral perturbation theory to the one-loop order. We also address the effects of fixing 
topology, by using saddle-point integration in the Fourier transform with respect to 6. 
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I. INTRODUCTION 



The low energy limit of quantum chromo dynamics (QCD) is full of nonperturbative 
phenomena such as quark confinement and chiral symmetry breaking. It has, however, been 
very difficult to analytically investigate these phenomena from the first principle due to its 
nonlinearity and strong coupling of interactions. Lattice QCD lj| and chiral perturbation 
theory (ChPT) have played prominent roles in studying such dynamical phenomena 

of QCD. The nonperturbative calculations of lattice QCD can be numerically performed 
utilizing the latest computational resources, while ChPT, an effective theory of pions, allows 
us to perturbatively treat the very low energy limit of QCD. They are complementary each 
other and have mutually developed. Lattice QCD in principle can determine the low energy 
constants of ChPT, some of which are difficult to determine from experimental inputs. On 
the other hand, ChPT provides a theoretical guideline of how to extrapolate the lattice data 
to the near-zero quark mass limit (chiral extrapolation) or the infinite volume limit (finite 
size scaling). 

An interesting extension of QCD is to introduce a CP violation term, known as the 9 
term. Since it is written as a total derivative in the QCD Lagrangian, it exists only when 
the gauge fields can have a nontrivial topological structure, or winding numbers. By partial 
integration, this 6 term can be regarded as the phase of superposition of different vacua in 
the Hamilton picture, which we call the 6 vacuum j^]. 

CP is invariant non only at ^ = but also at ^ = tt. It is, however, believed that CP is 

n 

spontaneously broken at = vr, where the theory has two CP violating vacua [5]. Moreover, 
it is expected that there exists a first-order phase transition in < < vr. These issues have 
been investigated mainly using effective theories 6|. 

It is little known in nature why this CP violation term is invisibly suppressed, which is the 
so-called strong CP problem. The neutron electric dipole moment has not been observed in 
the experiments, from which one can estimate ^ 10~^° 7j. The lattice QCD community 

remains, however, to be one of 
the most difficult problems in lattice QCD, since the 6 term gives a complex action, which 
leads to the notorious phase problem in Monte-Carlo simulations. With various elaborated 
approaches, many groups have investigated the 6 vacuum in the lattice simulations 

There exists another motivation in studying the 6 vacuum of QCD. One can consider a 
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fixed topological sector of the theory, by Fourier transform of the partition function with 
respect to 9. Expanding the vacuum energy as f{6) = C2^^/2 + c^O^/Al ■ ■ ■ and performing 
the saddle-point integration order by order, one can investigate how the topology affects the 
physical quantities and evaluate the difference between the physics between the 6 vacuum 
and the one at fixed topology 
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12l |. Inversely, it is also possible to determine the vacuum 



energy o: 



Refs. 
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QCD as a power series of 6, from the physical observables at fixed topology. In 



121 ]. a general formula to 0{1/V'^) is calculated which converts the observables at 
fixed topology to those in 6* = vacuum treating C2, C4, ■ ■ ■ as unknown parameters. It is, 
therefore, an important task to calculate these parameters within QCD or the low energy 
effective theory to quantify the effects of fixing topology. 

In fact, QCD at fixed topology can be investigated by lattice QCD simulations. Em- 
ploying the overlap Dirac operator 13|] which preserves the exact chiral symmetry, we are 



able to define the topological charge on^he lattice. JLQCD and TWQCD Collaborations 
are using a topology preserving way [14] to avoid discontinuities of the overlap fermion 
determinant, which considerably reduces the computational cost. The conversion formula 
between 6 vacuum and fixed topology is, therefore, essential in extracting the low energy 



constants such as the chira 



topological susceptibility [17]. 



condensate E 



15|, the pion decay constant F [l6| and the 



Because of its (//ofta/ topological nature, the 6 vacuum effect is totally infra red physics and, 
therefore, should be described by the lightest particles, or pions, within chiral perturbation 
theory. In this paper, we discuss ChPT with nonzero value of 6 as well as at a fixed topology. 
The formulation and general qualitative discussion are already given in Refs. js, Q. Our 
goal is to explicitly calculate the meson correlators to the next-to-leading order (NLO) in 
the 6 vacuum (and in a fixed topological sector) at finite volume, which may be directly 
compared with lattice QCD. In Sec. Hi] we will observe that the vacuum of chiral fields is 
located not at identity but at a nontrivial element on the SU{Nf) group manifold. In the 
p expansion of the chiral Lagrangian, this vacuum shift provides as a source which mixes 
the different CP eigenstates (Sec. IIII|) . In Sec. llVt we will calculate the 6 vacuum effect 
upto 0{e^) on the mesonic two-point correlation functions in ChPT to NLO. In Sec. IVl we 
will also address the physics at fixing topology, by using the saddle-point integration in the 
Fourier transform with respect to 6', as discussed above. The concluding remarks are given 
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in Sec. |VT1 



II. CHIRAL LAGRANGIAN TO THE LEADING ORDER 



We consider the A^j-flavor chiral Lagrangian in the 6 vacuum, 

C = ^i:i[d^U{x)^d^U{x)] - |Tr[7Wte-f?/iV/f/(a;) + U{x)^ e'^'^'f M] + ■ ■ ■ , (1) 

where U{x) G SU{Nf), S is the chiral condensate, and F denotes the pion decay constant 
both in the chiral limit. Here 6, the QCD vacuum angle, appears as the phase of the mass 
term, reflecting the picture that the 6 term can be converted to the chiral rotation of the 
quark bilinears through the anomalous Ward-Takahashi identity. In the mass matrix, 

M = diag( m^, m^/, • • - , 7711, m2, • • • ), (2) 
we have valence flavors and Nf dynamical flavors. 



In the partially quenched case [l9|], we use the so-called replica trick, where the calcula- 
tions are done with [Nf + + {N — A't,)]-flavor theory and then the replica limit N ^ 
is taken. The full theory results are precisely obtained by choosing = where irif 
denotes one of the dynamical quark masses. 

The system is assumed to be in the so-called p regime, where the Euclidean space-time 
volume V = L^T is large enough so that the perturbative expansion is performed according 
to the counting rule, 

d.r^Oip), a^)r^Oip), Mr^Oip'), T , L ^ O {I / p) , (3) 

in the units of a cutoff scale. 



Note that due to the vacuum angle 6', the expectation value of U {x) is located not at the 
identity but a nontrivial element of SU{Nf) (let us denote Uq). It is, therefore, useful to 
define the new variable and mass matrix, 

U{x) = UoU{x), Me = Ule'^l^fM, (4) 

and rewrite the Lagrangian, 

C = ^Ti[d,U{x)^d^iJ{x)] - ^i:i[MliJ{x) + iJ{x)^Me\ + ■ ■ ■ , (5) 
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where U (x) can be expanded around the identity as usual; 

where ^ is an element of SU{Nf) Lie algebra. 

By this vacuum shift, Eq.(l5]) explicitly shows that the nonzero 6 vacuum physics is equiv- 
alent to that in the 6 = vacuum but with a complex mass matrix Aig. Our task in 
the following is, therefore, to determine Uo (or equivalently Aie) and then to calculate the 
difference between the systems with the complex Aig and a simply real diagonal Ai. 



A. Vacuum shift Uq 



Let us first calculate the vacuum expectation value Uq, which minimizes the Lagrangian 
density of the zero- mode 



(7) 



For small 6, by parametrizing Uq = exp(?^°), the problem is equivalent to finding the 
minimum of a potential 



Nf 



Vie) = E 



sm 



'_e_ 

N 



o 



1 / e' 



+ ATr[e°],(8) 



where A denotes the Lagrange's multiplier to guarantee the traceless solution. 

It can be recursively shown, order by order, that all the off-diagonal parts of ^ cire zero 
For the diagonal elements, one obtains 

m 



a = 



-6 + a^r + Oie" 



where 



e 



m m ^ 1 



m 



Note that the partially quenched result is obtained simply by taking = rriy after the 
replica limit. It is also notable that = when m^'s are all degenerate. 
Now the original Lagrangian in Eq.Q is greatly simplified; 



6 \mi J rrij m'- ' 



(9) 



(10) 



Co 



Nf 



E COS ( — 6 — afO'^ 



(11) 
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from which one can read off the vacuum energy density at the tree level ^, as a function of 
6. In particular, the coefficients of 6"^ (the topological susceptibility Xt) and 6^ (we denote 
C4), 



Xt 



C4 



-m 



(12) 
(13) 



are important when we consider the effect of fixing topology as previously discussed in 



Refs. 
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12| 



Another important observation follows from the fact 

im{Me) = mei + o{e^) 



(14) 



where 1 denotes the Nf x Nf identity matrix. Noting Tr^ = 0, the contribution from the 
imaginary part of the mass matrix then becomes ~ 0{p^) so that we can neglect it at the 
leading order. ^ 

In this subsection, we have derived the vacuum expectation value Uq upto 0{6'^) level. 
The most part of this paper, however, requires only 0{6'^) contribution and one can neglect 
the 3rd term of Eq.Q or set = 0. 

B. Propagator of ^{x) 



Let us now expand the Lagrangian in ^, 



Here MfA9) is given by 



where mi{6) is defined by 



2tr 



m 



mi{e) = rrii cos —9 - aiO' ]=mAl + ^(^ 



mi 



1 m 



^2 



2mj 



(15) 



(16) 



(17) 



^ For the Nf = 2 case, a nonperturbative expression of the vacuum energy density is known [20'|, from 

which one can read off Xt, Ci [2l|, and any higher order coefficients C2ri's. 
^ It of course gives contributions at NLO. 
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and again there is no significant difference even if we extend the theory to the partially 
quenched one; we just set nii = rriy. 

The Feynmann propagator (in a finite volume) is then obtained: 

{^^A^)^kl{y))^ = SuS,,Aix - y, Mf^ie)) - 5,,5mG{x - y, M^^{9), M',,{9)), (18) 

where the second term comes from the traceless constraint Tr^ = 0. The definitions of A 
and G are given by (unless Nf = 0) 

1 pipx 

^(^■M^) - yZ-^, (19) 

1 pipx 

G{x, Ml Ml) = - E ^- ^, (20) 



where the summation is taken over the 4-momentum 

p = 2n{nt/T,n:^/L,ny/L,nz/L), (21) 

with integers n^'s. 

III. ONE-LOOP CORRECTIONS 

Since we are interested in at most one-loop corrections to the two-point functions, we can 
take, in advance, some Wick contractions in 0{p^) or 0{p^) contributions in the expansion 
of the leading order Lagrangian Eq. (fT5l) . 

Before performing the one-loop calculation, we here introduce the NLO terms of the 
chiral Lagrangian. Without source terms, we have 8 additional NLO terms whose low- 
energy constants are denoted by Lj's {i = l---8)'^. At 0{p^) and 0{p^), the terms with 



Li,L2,L3 (and Wess-Zunimo-Witten term 22|) do not contribute to our calculation. For 



the terms with L4, L5, Lg, the NLO correction is obtained in the same way as the 6 = case 
except for the change in the mass matrix; Ai KeAig, while the Lj and Lg terms require 
a special care of the imaginary part of the mass matrix, ImA^g = rhOl. 



^ When Nf = 2, due to the pseudoreahty, the number of independent terms is reduced to 5, of which 
coefficients are denoted by k's 
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Expanding the chiral field according to Eq.Q, we obtain 
1 



NLO 



E 

+E 



F2 



8Ly 

' F2 



■1 

_2 

Nf 



L,J2Mjf{e) + L,Mi^{e) 
f 

16 



X 



16L« 



where = 2mS/F2 and we have used [{ReMe)%i = [iMe=oY - ilraMe)%.i + 0{9'^). The 
first 4 terms can be absorbed into the redefinition of kinetic and mass terms, as usual in the 
case with ^ = 0. The 5th and 6th terms represent a peculiar contribution due to nonzero 6. 
We have kept the constant (but 9 dependent) terms in the last 3 terms for the calculation 
of topological susceptibility, which we will address later. 

A. Vacuum shift at one-loop 



Let us start with 0{p^) terms, which appears only in the nonzero 6 vacuum: 



= 9 



V2 



V2 



Nf 



Taking contractions of ^'s in the first term, and noting Y^iia = 0, it becomes 



where 



^ ml 
Bii = —X — 
rrii r ^ 



^ 9 



'Nf 



X: A(0, Mj) - G(0, Ml Ml - 16(iV^L7 + L,)Ml 



E I E ^(0, ^/) - G(0, M^, M^.) - lQ[NfLj + 



(23) 



(24) 



(25) 



which gives an 0{p'^) contribution. Here 9 dependence of the masses is dropped, since it 
gives only 0(6'^) x 0{p^) contributions; we have simply set Mf^ = Mfj{9 = 0) in (l25l) . Note 
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also that Ei MlBa = 0. 



The above hnear term in ^ in Eq. fl24p requires further shift in the vacuum Uq (xe *^/^/ ) ; 

[/,e-^^/^/ ^U'o^ diag(e-'^\ e-''\ ■■■), (26) 
where the phase of ith diagonal component is given by 

771 

0. = _0 - aiO^ -biO, (27) 
rrii 

with 

- ^/ 

b^^-B^^ + ^J:Bff. (28) 

Again, we have ignored 0{6^) x 0{p'^) contributions here. 
The meson mass is also shifted as 

M,5 (9) ^{M[.{e)f = ^ [nii cos 9., + m,- cos ) , (29) 
and the effective Lagrangian up to 0{p^) then reads 

+^M'9(Tre-3j:{a),^.}j. (30) 

B. Inserting sources 

Next, we consider insertions of the pseudoscalar and axial vector sources, p{x) and a^(x). 
Since the parity symmetry is broken by the 9 term, we will see that these source terms have 
unusual contributions which look like scalar or vector operators. It is therefore convenient 
to define the shifted Hermitian sources as 

p+ix) ^ i(f/;V(x)+p(x)f/^), (31) 

P'ix) = '-iU',^pix)-pixM), (32) 

o-t^^") = \iUo^(^^ix)Uo + a^.ix)), (33) 

= ^{Uoai,{x)Uo- a^{x)), (34) 



where we have assumed the original p{x) and a^(x) are both Hermitian and traceless matri- 
ces. In the following, we consider only charged meson type sources which have two different 
flavor indices. For this case the absence of the diagonal parts: [plu = [p~^]ii = [p~]ii = and 
= = = (fo^^ simplifies the calculation. 



C. One-loop effective Lagrangian with sources 



In the expansion of the leading Lagrangian Eq. fllSI) . we also have 0{p^) terms, 
1 



1 



(35) 



rom which contribution can be calculated in a straightforward way as in the case at 6 = 



Collecting all contributions so far, the LO + NLO effective Lagrangian, including the 
pseudoscalar and axial vector sources, is given by 



^ id 



Nf 



Nf 



+E 



Nf 



F3 



Nf 



+ ix)]ji[df,^{x)^{x) - ^{x)d^^{x)]i 



I Nf 



2F2 



X) 



A{0,Mm) 
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-E^^m/cos^/j - 4:L6{Y,Mjf{9)f + 4:Nf{NfL7 + L8)M^9^, (36) 

where we have omitted ^-independent constants. We also have omitted multi n point vertices 
for n > 3, which are irrelevant for the two-point correlation functions below. 
In the above result, the Z factors are given by 



1 



2F2 



E7^(A(0, M!,{0)) + A(0, Mf,{0))) ^ A{0, MUe), Mf,{9)) 



6 



+ 



Nf 



-8[L,Y.Mjf{e) + L,M^j{9) 



(37) 



z%{e) ^ 1 + 



2F2 



Nf 



G{0, Ml{e)) - 8(L4 - 2L,) ^ M}^{e) - 8{L, - 2Ls)M^.{e) 

f 



(38) 



Z'^{9) = 1 



2F2 



E;^(A(0,M^(^^)) + A(0,iW|/^))) 



+ A{0,MU9\MU9)) 



Nf 



L,Y.M}^{9)-rL,Ml{e) 



where 



A{x,MlMl) ^ G{x,MlMl) 



G{x,MlMf^+G{x,Mf.,Mf.) 



Ill jjj 



(39) 



(40) 



and its derivative 



dlA{x,MlMf^) = MlG{x,MlM\ 



3J' 



MlG{x, Ml M|) + M?.G(x, M2, M^) 



(41) 



are UV finite at a; = 0, and both vanish when = M|. 

On the other hand, A(0, M^) and G(0, , M|) are logarithmically divergent, and the 
divergent parts are evaluated by the dimensional regularization at D = 4 — 2e as 

A(0,M2) = 



- -h 1 - 7 -h In 47r 



167r2 Ve 



G(0,M2,M2) = 



167r2 



Ml + M| 



1 



Nf 



]^EA^//W| (7 + l-7 + ln47r)+---,(42) 



where 7 = 0.57721 • • • denotes Euler's constant. As is the usual case, these divergence can 
be removed by the renormalization of Lj's as 



miisub) - (^^ + 1 - 7 + ln47r - Iniil^^ , 



(43) 
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where L^^^subYs denote the renormahzed low energy constants at the subtraction scale fisub 
and 



1 Nf 1 1 1 \ 1 Nf 2 , 

74=g, 75 = ^, 76=g(2 + ]v^j' 77 = 0, 78=g(^^-^J. (44) 



As a result, Zp{9), Z]^j{9), and M-j{9) are kept finite, while Z^^^O) is still divergent but it 
does not affect the physical observables. 

After this procedure, one can replace A(0, M^) by 



A^(0,M^/iL) = j—,ln^ + g,{M'), (45) 



where gi denotes the finite volume contribution 



where Ki is the modified Bessel function and the summation is taken over the four-vector 
fl/i = n^iL^ with Li = L {i = 1,2,3) and L4 = T. Numerically, truncation at < 5 
already gives a good accuracy when ML > 3. For the exphcit expression of (^(0, , M|), 
see Appendix |Al 

It is noted here that the 8th term of Eq.f p6|) does not contribute to the NLO two-point 
functions, since we consider only off-diagonal sources, p{x)i^j or a^{x)i^j, whose coupling 
to the diagonal element is of higher order [next-to- next-to-leading order (NNLO)]. 



IV. TWO-POINT FUNCTIONS 



Pseudoscalar and axial vector correlators are obtained by the functional derivatives of 
the partition function 

Z\p,a^) = y'd^e-/'''"^'//(P''^^\ (47) 

with respect to the sources p and a^. 

Here we derive the two-point correlation functions of these operators in an irreducible 
representation which consist of two different valence quarks. We consider the most general 
nondegenerate case where their masses are denoted by m„ and m^,/, respectively. More ex- 
plicitly, we calculate the three types of correlation functions with zero momentum projection 
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(three-dimensional integral) , 



6 



+ 



6 



6 6 

+ 



6p{x)yy' Sp{x)y'yJ \6p{0)yy' 6p{0)y'y 



In Z\p, 



p=0,a^=0 



AoV{t, rriy, rriyi 
- fd'x 



6 6 
+ 



6 6 
+ 



6ao{x)yy' 5ao{x)y>y I \6p{x)yyf Sp{0)y 



In Z\p, 



(48) 



(49) 



p=0,a,j=0 



AoAo{t,my,my') = 



d'^x 



6 6 
+ 



5 6 
+ 



In Z\p,a^) 



(50) 



p=0,ap=0 



^5ao{x)yyt 6ao{x)y>v J \6ao{0)vv' 5ao{0)v'v, 
where we denote t = x^. As mentioned before,the partial quenching is performed by the 
replica trick: extending the number of flavors A'^^ — Nf + Ny + {N — Ny) then taking the 
limit ^ 0. 
Noting 



6 



5 



Sp{x)yyl 
6 

6ai,{x)yy' 



6p+{x) 



+ i 



6 



+ ^ ^ 



6a+{x), 



6P {X)yy' 

] — 



(51) 
(52) 



the correlation functions are given by 



VV(t,my,my,) = [Cpp{my,my,)] 



i-ioop cosh{M'J^°^{e){t-T/2)) 
M'J°°^{e) sinh(Mi-;^°°P(^)T/2) 



AVifr. r. )- (r. r. ^li-ioop ^^^MM.y-'-(g)(t - r/2)) 
AoVit,my,my,) - [CMm..my,)] 3inh(M-'-(^)r/2) 



(53) 



(gg-e) s^ 

4 F2 



(54) 



AoAo{t,my,my,) = [C%{my,my,)Y-'°°m'J°"^{e) 



cosh(M,V°°P(^)(t-r/2)) 



sinh(M„^;'°°P(^)T/2) 



) ^yyl /,N 



(55) 



Here, [Cjj;(m„, m^')]^~^°°P's denote the overall coefficients of which definitions are given in 
Appendix O 

Because of the CP violation, each correlator has a contribution from the two-pion state's 
propagation denoted by CJ^^t), C'-^' (t) and C^^'^^ respectively, which are of 0{e~'^^^*) at 
large time separation (See Appendix [C] for the details.). 
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The correction to the pseudo-Nambu-Goldstone boson mass is given by 



(56) 



and the pion decay constant is extracted from the PP(t, m„, m^/) and AV(t,my,myi) cor- 
relators in a standard way as 



= Fzf{e) 
1 



1 



d 



4F^ 9M2 



HUM' ) f^sub, 



M=M„ 



(57) 



The definitions of Hij{M'^) and H[j{M'^) are given in Appendix O 

One should note that for very small ^, the above formulas can be greatly simplified by 
Ignoring 0{e'^) X C(j92) corrections or just by setting W^O = and M^^,{e) = Mw(^). 



V. FIXED TOPOLOGY 



From the 6 dependence obtained so far, we can derive the correlators in a fixed topological 
sector of It is known that an observable at a fixed topology (let us denote Gq) is related 
to the one in the 6 vacuum [G'(^^)] by a formula 



i2i 



Gq = G{e = 0) + ^G{e)\ 



2xtV 



1 - 



C4 





XtV 2x'tV 



1 

+ ^G'(^)l^=o^ 



(58) 



which is valid in the general theories. Here, Xt = (Q^) denotes the topological suscepti- 
bility and C4 is the coefficient of 6'^ term of the vacuum energy of the theory. 

One should note that in ChPT, the 9 dependence only appears in the mass term, so that 
we can treat ~ C4 ~ O(A^) ~ OijP') in the p expansion. Therefore the factor l/^xtV) 
reduces the order of each contribution from fixing topology by 0{p'^), since 1/V = O(p^) in 
the p expansion ^. As a consequence, for nonvanishing G{6 = 0), one can easily calculate 
the "NNLO" contribution from fixing topology with the one-loop level calculation only. 



^ In the very vicinity of the chiral hmit (the e-regime), 1/xtV becomes 0(1) and cannot be treated as 
perturbation. Exact integrals over 9 is then needed, which is expressed by modified Bessel functions [18| . 
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A. Xt at one-loop 



Let us first calculate Xt within chiral perturbation theory at NLO 2l|]. For C4, as explained 
above, the tree level calculation we have given in Eq. flT^ is enough upto NNLO corrections. 
By switching off the source terms, one obtains from the mass term and terms in the last line 
of Lagrangian Eq. ( 136|) , 

1 



Xt 



1— loop 



^^^,.lnZ(p = 0,a, 



N 



rrii cos 9i 



0) 



1 - 



F2 



dd^ 



' Nf 




m m 
— + — 

rrii m 



Nf 



777 

A(0,Mj)-;^-G(0,M^,M^) 



■16L6 Ml - l6Nf{NfL-r + Ls)M^ 



(59) 



where all the masses in the last line are those at e = 0, or namely M^- = Mfj{e = 0). 
In the above calculation, we have used the fact that Si &i = [See Eq. (125]) ]. Note again 
that the UV divergence is precisely canceled by the renormalization of Lj's and therefore, 
one can replace A(0, M^) and Lj's with the renormalized values given in Eqs. fHSj) and fHSl) . 
respectively. Since the topological susceptibility is a coefficient of the 6^ term in the QCD 
vacuum energy, it does not depend on L4 and L5 at NLO, which only appear as corrections 
to the kinetic term. 



B. NLO correction from fixing topology 

As discussed above, the next-leading order correction from fixing topology can be calcu- 
lated at the tree-level. The above formula is then simplified: 

Gq = G{e = 0) + ^G(0)|,=o^ (1 - ^) + NNLO terms, (60) 

where Xt = X^'^ = ""^S is used. Furthermore, we can ignore all NLO x OiO"^) terms in the 
correlators G{9). Note here that we could have omitted term for small Q but kept it. 
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since it gives a (Q^) = xtV ~ 0{l/p'^) contribution when the topology is summed over again 
in the 6 vacuum. 

Substituting the expressions in the previous section into Eq. (pi]) , we obtain the correla- 
tors in a fixed topological sector of Q: 



{AoV{t,my,my>))Q = C^p{my,my>) 



M« sinh(M« T/2) ' 
sinh(MQ,(t-T/2)) 
sinh(M,^,,T/2) ' 
MZcosh{MZ{t-T/2)) 



smh{M^^,T/2) 



where the valence pion mass at fixed topology is given by 



m 



2xtV \mymy> 



XtV, 



and 



V Fzf{e = 0) ^ 



1 



m m 
— + 



1 - 



/ — 2 — 2 \ / 

/ m m \ 

+ 



XtV, 



AxtV \ml ml, 



XtV, 



-{FZl^ = 0))^ 



1 - 



m m 



1 - 



XtV, 



(61) 
(62) 
(63) 

(64) 

(65) 
(66) 
(67) 



AxtV Kniy rriy: , 

where Xt = 'mL. Note that each correction vanishes when summed over the topology, since 

m = XtV. 

From the above Q dependent correlators, the conventional extraction of the pion decay 
constant, of course, receives a correction from fixing the topology: 

C^p(mj,,m^/) 



w 



FZl^'{e = 0) 



1 



m m 



8xtV Krriy rriy: 
Note that the correction at NLO disappears when rriy = rriy'. 



XtV, 



(68) 



C. NNLO corrections from fixing topology 



In this subsection, we discuss NNLO corrections from fixing topology to the two-point 
correlators. Here we do not calculate two-loop diagrams at ^ = 0. They are already known in 



Nf =2 and 2-1-1 theories 



2j]-[3l|. Hereafter we denote them with a superscript "two-loop." 
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Ignoring the multipion states, the functional form of the correlators at two-loop in the 9 
vacuum has the following form: 



G{e,t) = C{9)f{M{9),t), 



(69) 



where C{9) denotes the time-independent coefficient, while f{M{9),t) represents the one- 
particle propagator with a mass M{9). 

In a fixed topological sector, however, the correction can not be factorized as x 
f{MQ,t) at two-loop level or more. Using notations 



1 



2X 



1— loop 

t 



V 



C4, 



Xt 



1— loop 



0'^"\9) 



gn 



-0(9), 



(70) 



for an arbitrary function O of 9, and directly substituting the above expression into Eq. (I58ll . 
one obtains at NNLO that 

92 



(71) 



M=M(e=0) 



where 



CQ ^ ci0) + C^'\0)6Q + C^'^^-f, 

^ M(2)(0)5q+ (^MW(0)+4M(2)(0)^i^^ 



(72) 
(73) 



Note that due to C*-^^(0)/C(0) dependence in , which is channel dependent, and the 
second derivative term, the correction in Eq. (17T]) cannot be simply absorbed into the mass 
shift in /(M,t). 



For more explicit expressions, we have to calculate the 2nd and 4th derivatives of various 
quantities with respect to 9. We summarize them in Appendix [Dl Using the notations there, 
the correlators at NNLO in a fixed topological sector Q are given by 



1 + [DUjJ'q^ + 2([^w]lo)^5^^ 



J J' 



(7,4) 

M=m2-/°°p(6I=0) 

vv' ^ ' 

[MfiWo^C + [{M^iVo + 4[Mi5],o[^7;]S) ?, (75) 
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where J and J' represent the operators P or A, 

cosh(M(t - T/2)) 



Msinh(MT/2) ' 
sinh(M(t - T/2)) 
sinh(MT/2) ' 
Mcosh(M(t-T/2)) 
sinh(MT/2) ' 



(76) 
(77) 
(78) 



and 



NNLO 



[C^=O(m,,m,0]'~^°°Px 



1 + ( [^pp]irio + 



0) ' Zr'(e = 0)/^ + f^^'^^^°'2 



+ 



(79) 



[C^p(m^,mv)] 



NNLO 



[Cy(m,,m,0]'"^°°Px 

+ [^ap\nlo 



zx.y'(^ = o) ' z):.-'(e = o)/ ^ apjlo 2 



+ 



[C^^(m^,m„/)] 



NNLO 



[C^t(m,,m,,)]'-^°°Px 

.'1(2) , 2[Z/](2) 



1 + [^aaIatlo + 



(80) 



(81) 



As is the mass correction, the decay constant at fixed topology is not uniquely ex- 
tracted from the correlators. If one has a good control of t dependence in Eg. (17^ on 
the lattice, however, a choice is to extract it from the coefficients [Cpp{ 

)]NNLO 



[Clp{'m^,m^,)] 



NNLO 



as 



NNLO 



[C^p(m^,m^/)] 



F[Z 



vv 



0)] 



2— loop 



1 + 



1 



2X 



1— loop 



1 - 



C4 



l(2) 



X < [Z^p]^^\o " '^I^p'pWlo + 



V \ XtV 2xtV, 



zf{e = 0) 



r yvv' ] (4) J- |- ,7^' 1 (4) 



PPjLO 



(82) 



In spite of a complicated channel dependence and nontrivial t dependence, one can see 
that the (axial) Ward-Takahashi identity is kept even at fixed topology; 



dt{AoV{t,my,m^>))Q = (m^ + m„/)(PP(t, m^/))( 



(83) 



which can be easily checked by starting from time derivative of the correlator in the 6 
vacuum, dt{AoV{t,my,mv>))Q, and performing integrals over 9. 
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VI. CONCLUSION 



We have discussed ChPT with a nonzero 9 term. As a result of CP violation, the vacuum 
of chiral fields is shifted to a nontrivial element on the SU{Nf) group manifold. We have 
calculated this vacuum shift at 0{9^) level, as well as the one-loop corrections, from which 
the topological susceptibility and C4, the coefficient of 9'^ in the QCD vacuum energy, are 
extracted. 

The CP violation also causes mixing among different CP eigenstates, between scalar and 
pseudoscalar, or vector and axialvector operators. We have calculated the mesonic two-point 
functions upto 0{9'^) to the one-loop order and 9 dependence of the pion mass and decay 
constant are obtained. 

We also have evaluated the effects of fixing topology, by Fourier transform with respect 
to 6. We found that the effect of fixing topology is considerably suppressed as expected; the 
tree level diagram only affects on the NLO corrections, one-loop diagram only contributes 
to the NNLO corrections, and so on. 

As applications of this study, it would be interesting to investigate three or four point 
functions, CP odd observables as well. It would be also important to compare our results 
with lattice QCD simulations. 

Acknowledgments 

The authors thank the member of JLQCD and TWQCD Collaborations for discussions 
and their encouragement to this study. In particular SA thanks S. Hashimoto and T. Onogi 
for discussions. HF thanks P. H. Damgaard for discussions and many useful comments. 
The work of SA is supported in part by Grants-in-Aid of the Ministry of Education, Cul- 
ture, Sports, Science and Technology of Japan (Grants No. 20340047, No. 20105001, No. 
20105003). The work of HF is supported by the Grant-in- Aid for Nagoya University Global 
COE Program, "Quest for Fundamental Principles in the Universe: from Particles to the 
Solar System and the Cosmos," from the Ministry of Education, Culture, Sports, Science 
and Technology of Japan. 



19 



Appendix A: Explicit expression for G{x, M^-, MJj) 

The diagonal part of the correlator 

G{x-y,Mm,Mf^{e)) (Al) 

can be, in principle, expressed in terms of A(x, M^). In this appendix, we consider the 
most general case with arbitrary number of flavors. The UV divergence of G at x = is 
also discussed. Furthermore, explicit examples for the degenerate theory and nondegenerate 



Nf = Ni + Ng flavor case will be given. The similar discussion can be found in Ref. 32 |. 
In the following, for simplicity, we omit the argument 6 in Mfj[d). Therefore, Mf^ means 
M^j{9), unless explicitly stated. 
Let us first define a function 

f{t) = —y^^, (A2) 

where k denotes the number of "different" quark masses, and > 1 is the degeneracy of 
the ith mass, which satisfies Y.i n-i = Nj. Here, we have ordered the masses Mf^ < M^j^^ ^j^-^ 
for all i. Noting that f{t) is a monotonically decreasing function, 

>--^i:^J^<0. (A3) 

and 

lim/(M| + e) = oo, lim/(M| - e) = -oo, (A4) 
fit) < 0, for t < Ml^ 

fit) > 0, for Ml < t, (A5) 

one can show that an equation /(t) =0 has k — 1 different solutions (we denote them by 
t = Mfj), each of them satisfying 

Ml<Ml<Ml,^^,, il<i<k-l). (A6) 



Therefore, /(— p^) can be written in a rational form: 

Y{]iP^ + M]^) 
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and G{x, Mf^, Mfj) can thus be expressed as 



Gix,MlM: 



Nf 



1 



4"^ A(x, M2^) + A(x, Mf.) + C(**)9m2 A(x, 



where the coefficients Ai*''^'s, eic. are given by the residues of 



ff) 



(A8) 



(A9) 



-t + Mi,){-t + Mfj) U'-'i-t + Mj 
[or -(-t + M2.)/2(t) for C^")], at each pole. Note that = when is equal to any of 
the physical masses. 

Next, we consider the UV divergence of G{x = 0, Mf^, Mjj). By expanding the denomi- 
nator of Eq.( l20i) in terms of masses, the UV-divergent part of G(0, Mf^, Mjj) can be written 

as 

1 



G(0,Mil,Mf, 



NfV 



1 



Nf 



where the logarithmic divergence of the last line is canceled by a renormalization of Lj's 
as seen in Sec. Illli Note that the quadratic divergence from term is absent in the 
dimensional regularization. 



Here, we give some explicit examples. For the fully degenerate case, i.e., equal masses 
M"^ for all sea flavor /, the above expression for G is greatly simplified. 



g(x,m2,m: 



Nf 



Nf 



33 



n' 



(All) 



(ML-M2)9m.A(x,M2)J {Ml 

For the Nj = Ni + Ns flavor theory, where we have A*"; quarks of mass mi and Ns quarks 
of mass m^, the equation f{t) = is easily solved and one obtains 



G{x,MlM^^ 



31' 



Afe')A (x, M2) + 5p)A(x, M|) + 5f )A(x, Ml)] {Ml ^ M: 
A(")A (a;, M^) + B^''^A{x, Ml) + C^") A(x, M^.)] (M^. 



(A12) 
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where 



{Ml - M^^){MJ^ - M2) ' ^ (M2 - Ml)(Mf, - M|) ' 
J.i^J) ^ {Ml-Ml){Ml-Ml) ^ N,N,{M?^ - Mjf 

{M^-Mf,){Mfj-Mi)' N]m-M^r 

Ml-M^^ 

where M^ = 2m{L/F^, M^ = 2msJ:/F^ and M^ = {N^Ml + NiMl)/Nf. 



Appendix B: Correlators of ^'s at finite volume 

In this appendix, we hst several useful formulas for the correlation functions of ^ fields 
or A(x, M^) at finite volume V — L^T. In particular, we consider the zero-mode projection 
or the three-dimensional spatial integrals. 

A useful identity is 

E ' llrl - k.M)e--(^-/.) + ,(_,M)e-(^-/^)l , (Bl) 



r (2p^' + M2 4Msinh(^) 



which holds for an arbitrary regular function g(j>). For example, by setting g( = 1, it is easy 
to obtain 



dx/X[x,M)-^^ sinh(MT/2) ' 



2. lsinh(M(t-T/2)) 



d'xdoA(x,M') = ^—^ -V^, (B3) 

" ^ ' ^ 2 sinh(Mr/2) ' ^ ' 



where we denote t — x^. 
Rather nontrivial ones are 



C^2{t,Ml,Ml) = j d^xA{x,Ml)A{x,Ml) 



T cosh(|g;|(t-r/2))cosh(|g^|(t-r/2)) 
^?=(i:;.,.3) 2ki|smh(|gi|r/2) 2|g^| sinh(|g^|r/2) ^ ^ 



CaAA(t, M2,M|) = j d^xdoA{x,M^)A{x,Ml 



T sinh(|g^|(t-T/2))cosh(|g^|(t-T/2)) 
^.-=(i:^...3) 2sinh(|gi|T/2) 2|g^|sinh(|g^|r/2)' ^ ' 



CaAdA{t,MlMl) = j d^xd^A{x,Ml)d^A{x,Ml), 
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T 



E 



smh(|g;|(t-T/2))sinh(|g^|(t-T/2)) 



^ .-=fef:^2,.3) 2sinh(|gi|r/2) 2sinh(|g^|r/2) ' 
CA92A(^,Ml^M|) = ld^xA{x,M^)d^,Aix,Ml), 

^ cosh(|gl|(t-T/2))|g^|cosh(|g^|(t-T/2)) 

V 



(B6) 



(B7) 



,^=(,u,.,,s) 2|gl|sinh(|gi|T/2) 2sinh(|g^|T/2) 

where \ql\ = sjq^ + Mf. As one expects, they are 0{e {A^i+M2)t^ large T, which describes 
the two pion state's propagation. 
In this paper, we also need 

CAA^(t, Mh Ml Ml) ^ I cfx I d^yA{x - y, M^)A{y, M|)A(|/, M|) 



1 



:,ipot 



X 



Jo 



+ dx 



In a;, M|,M| 



167r2 



(B8) 



and 



CsAAA(t, Ml Ml Ml) = Jd'xJ dSdoA{x - y, Ml)A{x - y, Ml)A{y, M 

+ 1 - 7 + In 47r^ 



1 ^ ipoe'^"* 



rp ^ y 



327r2 Ve 



+ 



/ — x) 

Jo 



lnF{pl,x,MlMl) + l 
167r2 



.(B9) 



Here F(pg, x, M|, M|) = x{l - x)pl + xM| + (1 - x)M|, and 



h''ix,po,F) = y: 



-tp a 



= -gcos(p.(.-l/2)a„)^ 



M=\fF 



(BIO) 



where the summation is taken over = Note that /i^(x,po; is exponentially small 

~ 0[e~^^). Since both of In F(pq, a;, M|, M|) and F) ai'e regular with respect to 

Po for nonzero masses and are symmetric under the flip of po ~^ ~Po, one obtains 

I ^.^ IW 1/ I / 'f \ \ \ ^ 

C^A<t,MlMlMl) 



1 cosli(.Ui(f - 772)) 



2Mi sinh {MiT/2) 



1 



In 1/ 



167r2 



+ Mm2,m2,m|) 

(Bll) 
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and 



C9AAA{t,MlMl,Mi 



where 



sinh(M3(t - T/2)) 
2sinh (M3T/2) 



i + 1 - 7 + ln47r) 
32^^2 



+ h' {Ml Ml Ml 



(B12) 



h{MlMlMl 



dx 



In F{-Ml X, Ml Ml + 1 



- cosh (^Mi (^x - ao^ 



167r2 
d ( M 



h' {Ml Ml Ml 



U7r2|a 



■i^i(M|a|) 



M=y/F(-M^,x,M^,Mi), 



(B13) 



c/x(l — x) 



\nF{-Mlx,MlMl + 1 



167r2 



^ cosh ( M3 ( X - ^"j ao^ ^ 



M 



\47r2|a 



■Ki(M|a|) 



M=y/F{-M^,x,Ml,M^)^ 



(B14) 



Note that the logarithmic divergences are canceled by the renormalization of Lj's (See Ap- 
pendix [C]). 

Appendix C: Details of two-point correlators 



In this appendix, we list several equations which are needed to calculate Eqs. (15311 . (15 
and ( 1551) . 

We first define 



Xij{x) 



Nf 



^ Y.A{x, Mj) A(x, Mff) - A(a;, M^.)G(x, Ml M^^ 
f 

-A{x, M')G{x, M' Ml) - 2A(x, M^. )G(x, Ml M?.), 



(CI) 



Nf 



Y.[d,A{x, Ml)A{x, Mfj) - doA{x, Mff)A{x, M^)] 
f 

-A{x, M^^)do[G{x, Ml Ml) - G{x, M^, mJ.)] 
+doA{x, M^.) [G(x, Ml Ml - G{x, Mj, Mj.)], 

Nf 

J2[d',A{x, Mfj)A{x, Ml) + d',A{x, Ml)A{x, M]j) - 2doA{x, Ml)doA{x, Mf^)] 
f 

+2d^,A{x, M^M{x, Ml Ml) + 2A(x, Mf.)d'oA{x, Ml M] 



(C2) 



3]' 



AdoA{x,Mf^)doA{x,MlMl). 



(C3) 
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In the above expressions, we have omitted the argument 6 in the mass unless its difference 
from Mfj{9 = 0) is important. Note that Yq*"' is always finite and Yq-^ = when rrii = rrij. 

Using the formulas in Appendix |Bl we can perform the zero momentum projection in a 
finite volume. It is straightforward to obtain 

cm = / d'xX^.ix), Clit) ^ I £xY^\x), C%Jt) ^ I d^xWi^iix). (C4) 

For example, Cx(t) is explicitly given by 

cm = J d^xXijix) 

f 

1 /k-1 \ 



N. 



Y: Af^C^t, Ml Ml) + i?(^^)CA.(t, M^, Ml) + C^^^^dM^C^t, Mg, M' 



f \ I 



\M=M,, 



( Y: Ai''^C^<t, Ml Ml) + B^"^C^.{t, Ml Ml) + C^"^dM^C^.(t, Mj, 



2 



\M=M, 



f \ I 
k-l 



E Af^^C^.it, Ml Ml) + Ml Ml) + Bf^C^.{t, Ml Ml 



r2 /\>2\ I jyii-i) 



2 T\j2\ , nfe')/ 



N. 



f \ I 



r2 i\j2 ■\ 



(C5) 



where k is the number of different sea quark masses, and the coefficients y4,i?,C"s and Ml 
are given in Appendix |Al 

Rather nontrivial integrals are 



£x y AA(x-i/,Mj.)X,,(|/) 

1 cosh(Mij(t-T/2)) 
~ 2Mij sinh {MijT/2) 



Nf, 



+ 1 - 7 + In 47r 



+ H,,{M, 



I] J 



(C6) 



d'x J d'yd'zA{x - y, MlX,,{y - z)A{z, M, 
1 d_\ j cosh(M(t - T/2)) 



'2MdM) \ 2M sinh (MT/2) 



2 Nj^ 



+ 1 - 7 + In 47r 



H.jiM' 



M=Mi. 



(C7) 



and 



73. / ^4„,v^.v. „,^A.„, .^2, _ l sinh(M,,(t-r/2)) 2^ 

2 sinh (M,,T/2) 



d'x d^yY^\x-y)A{y,Ml 
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where 



+ 24^^)h{M\Ml,Ml 
{B^ + 2Bf'^)h{M\ Ml, MfJ + {B^"^ + 2Bf)h{M\ MJ, mJ.) 

M^=Mj. 



(C9) 



TV, 



/ 



+ 



AT. 



/ 



(/^'(M^ M^, M^) - h'{M\ Ml M^)) 
-B^"\h'{M\ Ml Mf.) - /^'(M^ M,^,, M^. )) 



+C(**)5m.(/^'(M2, MJ, M,2) - h\M', Ml Ml))U=M„ 
-C^^^^d,Ah'{M\ Ml Ml - h'{M\ Ml MM 



(CIO) 



Note that 1/e divergence is canceled by the renormahzation of L-j and Lg. 
Finally we present the overall coefficients in Eqs.( l53l) . ( l54l) and ( l55l) : 



[C|,p(m^,m^/)] 



l-loop _ (^ot' ^(^))^ vdd' 

— ^ ) 



{0\ 



(Cll) 
(C12) 
(C13) 



where we have defined S^^/°°^(6') = S(Z]^'(^)Z™'(0))^ and dimensionless coefficients Z's 
are given by 

1 + cos(^^ + e^,) M^e^ ( d 



M=M, , 



z%{e) 



2p2 I dM'^^'"''^^ ^l^sub) 

„ + e,,) {?,2{NfL\{^l,u^) + Llifisub)) + Hl,\Ml„yi,^,))} , (C14) 



M^e 

cos^„ + cos^„/ M^e^ I d 



2F2 aM2' 



M=M„ 



4F2 



4< 



(C15) 
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1 + cos(6't, — O^i] 



d 



2 2F2 \ dM^ 



HUM' 



M=M, , 



F2 



where Ll^^fisub), LKfisub) and 



Nf 



Hvv'iMyy', Hsub) = Hyyi{Myy,) + —r — — In/i^yfe 



(C16) 



(C17) 



represent the renormahzed values at a reference scale ^sub- In these functions one can ignore 
the argument of the masses and just put M^^, = M^^,{9 = 0), since they only appear in 
0{e^) X terms. 



Appendix D: 9 derivatives 

In this appendix, we list the 6 derivatives of various quantities, which are needed to 
evaluate those in a fixed topological sector shown in Sec. |Vl Here, we evaluate them at 
NLO for the second derivatives while at LO for the 4th derivatives. 

For the mass matrix, we obtain 



Mi, (9) 





1 w? 


-Mij 


y2mimj J 




m m 
rrij 





— 2 

m 



rriimj 



(Dl) 



(D2) 



where the argument ^ = is omitted. 

For the one-loop propagator, or the chiral-log terms, we have 



A(2)(0,M2) ^ -Ml 



— 2 

m 



rriimj 



9Af2A(0,M2)|M=M., 



(D3) 



Using this, it is straightforward to calculate 



G(2)(0,M|,MJ.) = ^G(0,M2(^),M^.(^^)), 



(D4) 



A(2)(0,M?i,Mj.; 



G(2)(0,M^,m2.) - -(G(2)(0,M2,m2) + G(2)(0,M2,m2.)), (D5) 



and 



2F2 



Nf 



G(2)(0, Ml, Ml) + 8(L4 - 2U) M 



2 

// 



m 



27 



+8(L5 - 2U)Ml 



m 



(D6) 



2F2 



TV, 



+8 (^^4i:^//^+^5M^. 

It is then easy to obtain the second derivative of the mass at NLO: 



mirrij 



.(D7) 



[M, 



(2)i 

ij \NLO 



M- 



1— loop 



(^ = 0) 



m 



M2 



\2rn.jm,j 
2F2M2. 

For the overall coefficients of the correlators, we have 

1 / ^ ^ X 2 



0) 



(D8) 



PP\LO 



[■^PPiNLO 



m ^ m 
2 Vm„ m^/ 



Mlo +{ — + — ]ibv + K') - ^ M'^dM^HUM', fisub) 



m 



m 



rriv rUv' 



F2 



+ 



m m 
— + 



{32(iV^L;(/i,„,) + + i/;,,(M^^, 

Trio,/ / ^ ' 



[Z. 



1 I m 



™2 



— 2 

m 



AP\LO 



2 ml, ^ 



(D9) 

(DIO) 
(Dll) 



1(2) 



M2 



F2 



\2. 



AAlLO 



r7OT'i(2) 
L^AAJ iVLO 



m 



m 



2 Vm„ m„/ 



M2 



M2aM2i/,V(M2, 



m 



m 



Their 4th derivatives are given by 



r yvv' 1 (4) 

[■^ppIlo 



r yvv' 1 (4) 
V-^APlLO 



m ^ m 
2 Vm„ m,/ 



— + 



KAMI) 



rh 



+ 24(a^ + a^/) 



+ 24 



m 

Qy H Qy' 



(D12) 
(D13) 

(D14) 

(D15) 
(D16) 
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